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In the context of loop quantum cosmology, we parametrise the lattice refinement by 
a parameter, A, and the matter Hamiltonian by a parameter, 5. We then solve the 
Hamiltonian constraint for both a self-adjoint, and a non-self-adjoint Hamiltonian 
operator. Demanding that the solutions for the wave-functions obey certain physical 
■ restrictions, we impose constraints on the two-dimensional, (A, (5), parameter space, 

, thereby restricting the types of matter content that can be supported by a particular 

^P' lattice refinement model. 
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INTRODUCTION 



o 

cn 

U: 

I , Loop quantum gravity is a canonical quantisation of general relativity based on a Hamil- 

^1 tonian formulation with basic variables the connection, which basically carries information 
about curvature, and the triad, which encodes information about the spatial geometry. Re- 
CS) ' ducing the dynamical variables of the full theory to homogeneous and isotropic models, one 
^ ! gets loop quantum cosmology [1], which is not a field theory. 

QQ \ The fundamental variables of loop quantum cosmology are the holonomies of the SU(2) 

■ connection, A\., [i refers to the Lie algebra index and a is a spatial index with a and i 
^ . taking values 1,2,3) along a given edge, and the corresponding conjugate momentum, which 
^ i is the flux of the densitised triad, E^, through a two-surface. Assuming spatially flat, 
I homogeneous and isotropic models, the connection is given by a multiple of the basis one 
O ' forms, and the triad is obtained from the determinant of the fiducial fiat metric, ^Qab, which 
defines the volume, Vq, of the elementary cell, V. All integrations are performed over the 
^ i fiducial (elementary) cell, V. 

I To proceed with the quantisation procedure, one has first to construct the Hamiltonian 

^ operator. Dynamics are then determined by the Hamiltonian constraint. We emphasise 

that while in the full theory there are an infinite number of constraints, in the reduced 
homogeneous and isotropic case there is only one integrated Hamiltonian constraint. Matter 
is introduced by adding the actions of matter components to the gravitational action. Thus, 
one just adds the matter contribution to the Hamiltonian constraint. One then obtains 
difference equations, analogous to the differential Wheller-DeWitt equations. 

Phenomenological reasons 0, 0, 0] require the parameter appearing in the regularisation 
of the Hamiltonian constraint not to be constant. Considering an underlying lattice which is 
being refined during dynamical changes of the volume, one allows the number of vertices on 
the closed loop making up the holonomies to vary dynamically. One has then to implement 
this requirement in the quantisation procedure of the Hamiltonian constraint. 

We parametrise the lattice refinement and the matter Hamiltonian by introducing two 
parameters, A and 5, respectively. Considering the Hamiltonian operator in the self-adjoint 
and the non-self-adjoint cases, we solve the constraint equation for both fixed and varying 
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lattices. Demanding that the solutions should satisfy certain physical requirements, we 
impose constraints on the two-dimensional, {A,S), parameter space. 



II. CONSTANT LATTICE 



The gravitational part of the Hamiltonian operator, Cgrav, can be written in terms of 
SU(2) holonomies, hi, and the triad component jj, determining in the fiat {k = 0) model 
the physical volume of the fiducial cell, V, as [2I, ^ 



2i 



grav 



h 



(mo)-1 



V 



sgn(]3) 



(2.1) 



ijk 



where k = 8nG, and V = Ip^/^ denotes the volume operator. Note that we use the 
irreducible representation, J = 1/2, since in this case the Hamiltonian constraint is free 
of ill-behaving spurious solutions [21, 0]. The holonomy, hl^°\ along the edge parallel to the 
i^^ basis vector, of length fioVQ^^ with respect to the fiducial metric, is [sj 



M^''' = cos 



(2.2) 



where 1 is the identity 2x2 matrix and = — i(Tj/2 is a basis in the Lie algebra SU(2) 
satisfying the relation 

r,r, = (l/2)e,jfcr'= - (l/4)(5,j , 

with cTi the Pauli matrices. The pair (c, p) denotes the coordinates of the two-dimensional 
gravitational phase-space. The triad component p determines the physical volume of the 
fiducial cell and the connection component determines the rate of change of the physical 
edge length of the fiducial cell. They are related through 



{c,p} = Y 



(2.3) 



with 7 the Barbero-Immirzi parameter representing a quantum ambiguity parameter of the 
theory. 

The action of the operator exp[i{fiQc/2)] on the basis states, with 

p\fi) = {K^h\fi\/6) \fi) , 

where /z (a real number) stands for the eigenstates of p, satisfying the orthonormality relation 

(/il|/i2) = 5/.Ji,/i2 ) 

reads 



exp 



l/i) = exp 



djj, 



1/^) = l/^ + /^o) ; 



(2.4) 



jiQ is any real number. 

The action of the holonomies, h\^''\ of the gravitational connection, on the basis states 
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is given by [5j 

/iS'^^V) = (csl + 2snri) l/i) , (2.5) 



where, 



Thus, 



cs|/i) = cos(/ioc/2)|/i) = [ l/i + /io) + l/i - fio) ] /2 , 

sn|/i) = sin(/ioc/2)|;u) = [ \fi + fio) - \fi - fio) ] / {2i) . (2.6) 



;^(Mo)^(Mo)^(Mo)-l^(Mo)-l|^^ 

= [(cs^ - sh^) 1 + 2 (1 - ATjTi) a^sn^ + 4 (r, - r^) Ira sn^] , (2.7) 



and 



h 



(mo) 



= (V - csVcs - snV^sn) + 2ri (^ral/sn - sn\>cs) . (2.8) 
Substituting Eqs. (EZD and (ESD into Eq. (O) we obtain 



CgravlAt) = 2fc^3 s ^^^h^ TsliV^cs - csysh") . (2.9) 

Using Eq. (12. 6p we recover the known expression for the action of the gravitational part of 
the Hamiltonian constraint, namely^ 

^ / ^ N 1/2 

Cgravl/^) = ^ l-^j f^o'Sif^) [1/^ + 4;Uo) " 2|/i) + |/i - 4/io)] , (2.10) 
where S{fi) is defined by 

S'(/i) = |/i + /io|^^^ - |/i-/Uo|^^^ • (2.11) 
To make the Hamiltonian operator self-adjoint we simply define 



^grav 2 



2 (pgrav ~\~ C^graJ^ ; (2-12) 



-'^ Being interested in the large scale behaviour of the loop quantum cosmology equations, we neglect the 
sign ambiguity that arises from the two different orientations of the triad. 
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which acts on the basis states as 



1 / h 



8 \6k'J'- 



1/2 



S{fl) + S{fi + 4/io) I + 4yUo) 



(2.13) 



Taking the continuum hmit {fi ^ /io) of the Hamiltonian constraint equation 



(2.14) 



and expanding the general state in the kinematical Hilbert space in terms of the basis 
states, l/i), as 

\^) = J2'^,{<P)\fi) , (2.15) 

where the coefficients are not continuous with respect to /i and the dependence of the 
coefficients on represents the matter degrees of freedom, we get 



n 



" 8 \6k^^ 
2/io 



1/2 



2fl, 



12/i 



2 (^M+4mo - 2^^. + + (/^o) 



(2.16) 



where Tit/il^I/) = '^^</)|^) is assumed to act diagonally on the basis states We note that 
the kinematical inner product of the general states reads 



(2.17) 



with the requirement that a state in the kinematical Hilbert space must have a finite kine- 
matical norm. The basis states are eigenstates of the volume operator and while the 
eigenvalues n are valued on the whole real line, the states are normalisable with respect to 
the kinematical inner product 

Assuming that the wave-function does not vary much on scales smaller than 4/io (known 
as pre-dassicality [§[), one can approximate \E'/i(0) as '^fj.{4>) ~ \E'(/i, 0). Then using Taylor 
expansion the constraint equation reduces to the Wheeler-DeWitt equation ^ 



7i:^$(/i,0) = 6 



6kY 



1/2 



(9/i^ 



+ /i 



(2.18) 



we have re-introduced the dependence on the matter degrees of freedom 0. The quantum 
constraint, which is a difference rather than a differential equation, constraints the coeffi- 
cients \E'(/t, 0) to ensure that |\E') is a physical state. 
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The non-self-adjoint version of the operator produces a different factor ordering, namely 



= 12 



h 



1/2 



(//2vl>(;,,0))+0(/io) 



(2.19) 



which affects the conditions on normalisability to be discussed later. 



III. LATTICE REFINEMENT 

The case of a dynamically altering holonomy length scale, is required for several 

phenomenological reasons [l,]!, 0]- However, this is not just a naive substitution, fiQ 
in the previous equations. One can immediately realise that this would lead to difficulties, 
since there would be extra terms arising in Eq. 02.181] as a result of the dynamics of the 
underlying grid. 

To derive the correct constraint ecmation we need to introduce the varying length scale 
into the definition of the holonomies |5| 



hi = exp 



—la-i _ 
——lie 



(3.1) 



where the reader should keep in mind that depends on /x. Geometric considerations [5| 
imply that, after quantising. 



exp 



/iC 



|^(/i,0)) = exp 



d_ 



(3.2) 



This however is no longer a simple shift operator, since /x is a function of fi. Consider 
changing the representation from /j, to 



(3.3) 



where /iq is a constant. In this representation we have 



hilu) = exp 



|z/) = exp 



dz/ 



\u) = |z/ + /io) • 



(3.4) 



We can then proceed as before and define 

sh|z/) = sin 
cs|z/) = cos 



/iC 

2 _ 


' ' 2i 




fic 
2 


' ' 21 


\u + flo) + |z/-/io) 



There is however a problem in defining the volume eigenvalue, since this requires an explicit 



6 



relation between v and /i given by /i. Assuming 



/i — /Xq/^ 1 



one has 



l-A 



^o(l-A) 

(up to a constant that can be set equal to 0), leading to 



(3.5) 
(3.6) 



V\v) 
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3/2/(l-A) 



(3.7) 



A. Non-self-adjoint case 
Let us calculate the action of Eq. fl2.ip on the basis state 



where 

and S{u) is defined by 



a = yUo(l - A)/fio 





3/2/(l-A) 




3/2/(l-A) 


(z/ + /io) a 




(z/ - /io) a 





(3.8) 
(3.9) 

(3.10) 



One can easily check that A = reproduces Eq. ( 12.10p . if it is taken that /io = /io- After 
a long but straight forward expansion in the z/ ^ /io limit and under the assumption of 
pre-classicality, one finds 



1/2 



^ -3/2(l-A) ^(1-4A)/2(1-A) 



1 - AA d^{iy, 0) (1 + 2A) {AA - 1) 



(z/, 



9z/2 (l-^)i^ du Ail-Afu"^ 
+0(/io) . (3.11) 

For a non-self-adjoint Hamiltonian operator, the Hamiltonian constraint equation reads 

(3.12) 



V dv 
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where 



l-AA 

(1 + 2A) {AA - 1) 

A{l-Af 
^3/2/(l-A) /6^^3y/2 



B 
C 

^ 12(1 -A)2 V h 

We note that for a fixed lattice, A,B,C and /? are given by 

/l = , 5 = 1 , C = -l/4 , (3={QK^^/hY'\iiolhf^ 112. 



(3.13) 



(3.14) 



Considering lattice refinement in the case of a non-self-adjoint Hamiltonian operator, one 
has u = fiQfi/fiQ. Thus, by keeping /io general all we have done is to re-scale /i. Setting 
/io = /io we get back Eq. fl2.19p . 

The specific lattice refinement yl = — 1/2[3|, Isjis clearly a particularly fortitious choice as 
it results in C = 0. Notice however that choosing A = 1/4 results in a further simplification, 
leading to 



+ 



12 



3\ 1/2 



(3.15) 



assuming fiQ = /io. This unphysical lattice refinement choice results in dynamics that are well 
approximated by the Wheeler-DeWitt equation in the large scale limit (for slowly varying 
wave-functions). This remark highlights the importance of understanding the origin of the 
lattice refinement in the full theory. Unfortunately, at present there is little theoretical 
reason for discounting such physically unacceptable scenarios, and one must rely on his/her 
phenomenological intuition. 



B. Self-adjoint case 

Let us repeat the above procedure for the case of a self-adjoint Hamiltonian operator, 

"^grav — (^^grav + ^grav)/^ ' 

Acting on the state one has 

1 / h 



8/iQ \6k7 



S{u + Aj2o) + S{u) \\iy + 4j2o) 



-AS{u) + |^(z/ - 4/io) + S{u)^\u - 4/io) 



Expanding Tigravl^) = —'H^\'^), one obtains 



92^(z/, 



+ 



B (9^(z/, 
2^ 



(3.16) 



(7z/-2^(z/, 0) + /?7^<^z/-^/2^(z/, 0) + O (/io) = , (3.17) 
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where B and /3 are given by Eq. (13. 13b ) and Eq. (13.13b ). respectively, and 

1 - lOA 



B 



1-A 



^ ^ (1 + 2A)(4A-1) + 12A(2A-1) 

8(1 -A)2 ■ ^ ■ ^ 

As expected, setting A = and /io = /io gives back Eq. (I2.18p . Once again we see that the 
choice A = —1/2 produces a particular simplification. 



C. Physical sector 

In general, not all solutions to the quantum constraint equation, Eq. (I3.17P in the case 
of a self-adioint Hamiltonian operator, are normalisable with respect to the physical inner 
product [5|,|y|. In what follows, we are only interested in physical states. The physical Hilbert 
space consists of solutions to the quantum constraint equation which have finite norm with 
respect to the physical inner product. The inner product on physical states can be obtained 
by requiring that real classical observables be represented on the physical Hilbert space by 
self-adjoint operators [6|. The physical inner product has been calculated js], @] if the only 
matter source is a massless scalar field. Following the same procedure, we will compute the 
inner product of physical states for the model we are considering here. 

The (total) Wheeler-DeWitt constraint equation reads 

(^grav + 7i:<^) ^ = . (3.19) 

Since we are interested in the large scale limit, we approximate the matter Hamiltonian, H.^, 
with = z>^e (0) (the reader is referred to the next Section). Thus, 

e (0) $ = e (0) ^ = -zz-'TTfgrav^ • (3.20) 

In the classical theory, e (0) is a Dirac observable since it is a constant of motion Even 
though 1/(0) is not a constant of motion, assuming that z/(0) is a monotomic function (with 
respect to 0), then z/(0o) is a Dirac observable for any fixed 0o [sl]. Modulo an overall scaling, 
the unique inner product which makes these operators self-adjoint is |5| 

(^l|^2)phys= / du\u\'^,^2 . (3.21) 

The finite norm of the physical wave- functions, defined by Eq. (13.2ip . is conserved, i.e., 
independent of the choice of = 0o. From Eq. (I3.2ip one concludes that the solutions of the 
constraint are normalisable provided they decay, on large scales, faster than One 
arrives to the same conclusion for the case of a constant lattice, with u replaced by fi. 

It is important to note that, in general the approximation of Ti,^ = u^e (0) is only valid 
on the large scale i> limit, implying that the integrand of Eq. fl3.2ip is only valid for ^ 1. 
However, it is certainly necessary that the large scale behaviour of the wave-functions be 
normalisable with respect to Eq. fl3.2ip . Thus, the constraint we have found is a necessary 
but not sufficient condition for the wave-functions to be considered physical. 
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IV. SOLVING THE CONSTRAINT EQUATION 



To solve the constraint equation one needs to know the specific form of 7i^. In general, Ti^ 
has two terms with different scale dependence, however since we are concerned only with the 
large scale limit, one of these terms will be the dominant one. Making this approximation, 
one can write 

= e^{<j))fi'^ or = e^v'^ , (4.1) 

where the functions e^, e^, are constant with respect to /i, z/, respectively. The general 
analytical solutions read 



non— self —adjoint, fixed lattice S^^y^ 



self— adjoint, fixed lattice 



C\ ^2/(3+25^ 
+C2 1^2/(3+25^) 

'-'1/^ -'^3/(3+25^) 



3 + 25, 



_Vf^^(3+25^)/4 



3 + 2(5, 



4./e 



3 + 25 



^(3+25^ )/4 



non— self —adjoint, varying lattice v]/(jy) 



self —adjoint, varying lattice v]/(jy) 



+C2/i^/^ ^^3/(3+25,) 
C<^^^3A/2/(A-l) J^^^^_^ 

^^^(l+8A)/4/(l-A) J^_^^ 



^ ,,(3+25^)/4 



(4.2) 



with 



25,(1 -A) + 3 
4(1 - A) 

v/3(12A + l) 
4(1 -A) ' 



(4.3) 



J and y are Bessel functions of the first and second kind, respectively, and Ci,C2 are 
integration constants. Note that we suppressed the dependence for clarity. We wrote 
explicitly the solutions for the non-self-adjoint, as well the self-adjoint case for both a fixed 
and a varying lattice. In particular, for the p/iT/szca//?/ justified choice A = —1/2 [3I, 0] the 
solution of the Hamiltonian constraint equation, in the case of a non-self-adjoint Hamiltonian 
operator and a varying lattice, reads 



2^/e 



IJ+Tj V 5,^ + 1 



^^(5.+l)/2 



(4.4) 
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Among such solutions, we will only consider the physical ones. This immediately eliminates 
the solutions to the non-self-adjoint Hamiltonian constraint, since it is not possible to find 
self-adjoint Dirac observables for these cases. We nevertheless find interesting to compare 
the non-self-adjoint solutions to the self-adjoint ones and we thus apply the norm defined by 
Eq. fl3.2ip to both (self-adjoint and non-self-adjoint) sets of solutions. This is done simply 
to complete the formal comparison between the two cases and it is to be remembered that 
the normalisation constraint produced here is only rigorous for the self-adjoint case. 
Using the asymptotic expansions of the Bessel functions, 

r T ( ^ [Y f Pn 71 
hm jRiz) \ — cos 2; 

lim Yg(z) \ — sin ( z — — 

z^oo ^ \ TVZ \ 2 4 

we find that the solutions oscillate within an envelope that scales as 

non— self —adjoint, fixed lattice v]/(^) qq ^^{3+25/j)/8 

self —adjoint, fixed lattice \[/^^^ ^~{1+2<5^)/8 

non-self-adjoint, varying lattice v]>^^^ q(_ ^[12A-25^(1-A)-3]/8/(1-A) 

self-adjoint, varying lattice ^ ^ ^[16A-l-25i,(l-A)]/8/(l-A) _ 

As shown in Section [111 CI if the solutions are to be normalisable, ^'(z^) must not grow faster 
than \E'(z/) oc u~^^^'^^\ which imposes constraints on the scale dependence of the allowed 
matter component. More precisely, 

1 

non- self- adjoint, fixed lattice 

self-adjoint, fixed lattice 

non- self- adjoint, varying lattice 

self-adjoint, varying lattice . (4.6) 

In the self-adjoint and non-self-adjoint Hamiltonian operator cases with lattice refinement, 
the growth is taken w.r.t. u. If we require the semi-classical wave-functions not to grow w.r.t. 
fi, there is an additional constraint, namely A < 1, to ensure that increasing fi corresponds 
to increasing i^. 

One should also keep in mind that since the solutions are only valid on large scales, one 
must ensure that the large argument expansions of the Bessel functions apply in this limit. 
The expansions are valid for 

5m > -3/2 , (4.7) 

in both cases of a self-adjoint and a non-self-adjoint Hamiltonian operator considering a 
fixed lattice, and for 

> . (4.8) 
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Non-self- ad joint, varyine; lattice , , , . 



g 1 a-VT? 3 



2 8 2 



Self- adjoint, varying lattice 

e((p) = 



/ // / /// // / / / Normalisable 

^' xW^ "^.^ ,\^?" W Non-normalisable 

■->. Normalisable only for Ci— 0orC2— 



FIG. 1: Regions in which the solutions to the different vacuum Hamiltonian constraint equations 
are normalisable at large scales. Note that seif-adjomt, varying lattice ^(^j^) oscillatory only for 
A < -1/12. 



in both cases of a self-adjoint and a non-self- adjoint Hamiltonian operator considering a 
varying lattice. Beyond these limits the wave-functions, on large scales decay like 1/y, 
where the argument of the Bessel functions is n^, or u^, respectively. Whilst these wave- 
function may be normalisable, they lack a semi-classical interpretation and hence would 



not produce classical cosmology at large scales [10|. The regions where the different wave- 
function coefficients are bounded are shown in Fig. [2] along with the limit of the expansion. 

A particularly interesting case is that of the vacuum, where Ti,^ = 0. This corresponds 
to e (0) = or 5,^ = —oo, which makes the norm calculated in Section [111 CI trivial. In this 
case the norm can be taken to be, 



(^l|^2) = / , (4.9) 

J 4>=(f>o 

and correspondingly for fi. The solutions to the four constraint equations are 

noil— self —adjoint , fixed lattice \[/^^^ Ci^^^'^ -\- C2/^~^^^ 

self— adjoint , fixed lattice \[/^^^ Cifl^^^^'^^^ + C2fJ'^^~^^^^ 

non— self— adjoint , varying lattice v|/^j/) ^{4^~l)/2 _|_ ^ ^(2A+l)/2 

self-adjoint , varying lattice v|>(^i/) = ^(§-4+1+ V36A+3) /4 _|_ ^^^(8A+1- V36A+3)/4 j^g^ 

Clearly the two cases of Hamiltonian operator for a fixed lattice are bounded only for specific 
choices of the integration constants, which amounts to special initial conditions. For the lat- 
tice refinement case however, there are several regions in which the solutions are bounded, 
shown in Fig. [H It is also worth noticing that only the self-adjoint lattice refinement case 
produces oscillatory solutions (for A < —1/12) and hence have a simple semi-classical dy- 
namical interpretation lo| . 
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V. LARGE SCALE CLASSICAL BREAKDOWN 

The form of the wave-functions indicates that the period of oscillations can decrease as 
the scale increases, which implies that at sufficiently large scales the assumption that the 
wave-functions are pre-classical may break down [l^, llj . This would then lead to quantum 



gravity corrections at large scale (classical) physics. The need to avoid this undesired event, 
was indeed one of the motivations behind lattice refinement. However, the level of lattice 
refinement necessary depends on the matter content. 

One way to investigate this is to look at the separation between the zeros of the wave- 
functions [4] (for an alternative method see Ref. [3|). For the non-self adjoint Hamiltonian 
case we find that the n*'* zero (for large scales) occurs at 

Tr(9^(A q\-|4(A-l)/[254A-l)-3] 

If A ^ in + Cf''-'^/^''^^^-''^-'^ , (5.1) 

4(A — l)^/e^ J 

where 

1 -1/^ c'A , (-"-1) ,11 



is a constant. Since Eq. (15. ip is derived from the large argument expansion of Eq. (14. 2p . it 
is only valid for 6^, > 3/2 /{A — 1). Using a Taylor expansion we find 

77 (4-2ii.)(l-A)-3 / 4(l-,5^)(l-A)-6 \ 

lim Aun = —^v ^^^^ +0[v 4(1-4) . (5.3) 

largei/ V / 

Note that the Taylor expansion is valid for 5 > 1 + 3/2/(^4 — 1). Using 



/io(l-^) ' 
we find 

lim Az/„ OC ^(f-2)(A-l)/2-3/4 _^ ^ / (5.-l)(A-l)-3/2N _ ^5 4) 

largei/ 

The lattice refinement will support all oscillations of the wave-function, provided Az/„ is 
larger than z/g, the scale at which the underlying discreteness becomes important. The 
condition for the continuum limit to be valid, is that the wave-function must vary slowly on 
scales of the order of [3| 

/ic = 4/i , (5.5) 

or, equivalently, 

fic = 4/io/^^ • (5.6) 

Then from Eq. (13. 6p one gets 

Equations (15. 4p and (15. 7p imply that lattice refinement will be sufficient to prevent quantum 
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All constraints 
satisfied 



FIG. 2: The regions of parameter space in which the wave-functions of the self-adjoint Hamilto- 
nian constraint equation with lattice refinement are physically acceptable. Notice that there are 
regions (crosses) in which the Taylor expansions used to calculate the large scale behaviour of the 
wave-functions are no longer valid, and hence whilst we can say that these wave-functions decay 
sufficiently quickly on large scales to be normalisable and are physical (i.e. oscillating), we cannot 
be sure that there is no new large scale behaviour due to the underlying discreteness. 



corrections becoming significant at large scales provided 

/(A^.)-A^+(|-2)A + i-|>0, (5.8) 

with further restrictions on ejy(0) for the case of equality j^. This is shown in Fig. [2J 
A similar calculation for the self-adjoint case gives 

Vr (4-2a^)(l-A)-3 / 4(l-g^)(l-A)-6 \ 

lim Ai/„ = u 4(i-A) +0 iu 4(i-A) \ ^ (5.9) 

large i/ ^ 

which is precisely the same equation we had for the non-self-adjoint case, Eq (15.31) . In 
addition, the Taylor and Bessel expansions used are valid in the same ranges as those of the 
non-self-adjoint case. This is not surprising since making the constraint equation self-adjoint 
is inherently a quantum operation and their classical limits should be the same. Thus, Fig. [2] 
applies to self-adjoint, as well as to non-self-adjoint Hamiltonian operators, in the lattice 
refinement case (albeit with a different constraint coming from the requirement that the 
coefficients be normalisable). 

By considering the underlying origins of lattice refinement, we can further restrict the 
allowed range to0<y4<— 1/2[3|. This allows us to examine the types of matter that 
cannot be supported by a particular lattice refinement model. The relevant section of Fig. [2] 
is replotted in Fig. [31 Notice that Eq. (13. 6p is needed to find the scaling behaviour of a 
particular matter component with respect to the scale factor, i.e. matter scaling like u^" , 
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FIG. 3: The full loop quantum gravity theory allows only the range < A < —1/2. Within 
this range we can see that the acceptable types of matter content are significantly restricted. In 
addition, notice that for a varying lattice {A ^ 0) it is not always possible to treat the large scale 
behaviour of the wave-functions perturbatively (dashed line with crosses). 

scales with respect to the scale factor as a'^^''^^~^K 

VI. CONCLUSIONS 

We have derived, in the continuum limit, the Hamiltonian constraint of loop quantum 
cosmology for a general lattice refinement scheme of the form p, = fJ,ofJ^^, for both the self- 
adjoint and non- self-adjoint Hamiltonian operator cases. We solved the resulting Wheeler- 
deWitt like equations and discussed the requirements the solutions must satisfy in order 
to be physically viable. These requirements give us constraints on the type of matter that 
can be supported by a particular lattice refinement model. We considered the following 
three requirements for the wave-functions: (i) that the coefficients of their basis expansion 
be normalisable, (ii) that they have oscillating large scale solutions so as to ensure that 
classical dynamics can be recovered and (iii) that they are pre-classical at large scales. 
Combining these conditions significantly constrains the allowable region of parameter space. 

In particular, for the case of a constant lattice, physical wave-functions are produced only 
if Ti.^ scales faster than and slower than a; an extremely severe restriction, given that 
many types of matter scale beyond this range. In the most popular lattice refinement model, 
A = —1/2, this range is extended so that physical wave-functions are produced provided Ti^ 
scales faster than and slower than a^, although it is not possible to treat the large scale 
oscillations perturbatively over a third of this range {a~^ a°). 

As a concrete example, an inflationary scalar field (i.e., one in which the potential term 
dominates over the kinetic term in the matter Hamiltonian), scales like a^. From our general 
procedure it is clear that this has a large scale breakdown of pre-classicality for the fixed 
lattice case, whilst this problem is resolved for the common, A = —1/2, lattice refinement 
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case, as was shown in [J] . This provides a further demonstration of the importance of mod- 
elhng lattice refinement in loop quantum cosmology, if physical results are to be produced, 
and it does so for a large class of such models. 

It is important to note that lattice refinement could, in principle, be much more compli- 
cated than the power law form {fl = yUo/^"^) used here, however even with this simplifying 
assumption the qualitative behaviour of different lattice refinement models is clear. In par- 
ticular, we have shown that the continuum limit of the Hamiltonian constraint equation is 
sensitive to the choice of model and that only a limited range of matter components can be 
supported within a particular choice. This further emphasises the need to support effective, 
phenomenological lattice refinement models with a deeper understanding of the fundamental 
theory. 
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